Ferromagnetic and random spin ordering in diluted magnetic semiconductors 
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In a diluted magnetic semiconductor system, the exchange interaction between magnetic im- 
purities has two independent components: direct antiferromagnetic interaction and ferromagnetic 
interaction mediated by charge carriers. Depending on the system parameters, the ground state of 
the system may be ordered either ferromagneticaUy or randomly. In this paper we use percolation 
theory to find the ferromagnetic transition temperature and the location of the quantum critical 
point separating the ferromagnetic phase and a valence bond glass phase. 
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I. INTRODUCTION 



In diluted magnetic semiconductors (DMS), e.g. 
Gai_2;Mn3.As, In1_3.Mna.As etc, long-range ferromag- 
netic ordering of the impurity (i.e. Mn) local moments 
is induced by the carriers (holes in most cases) which 
are contributed by the dopant impurities serving the 
dual role of acceptors and magnetic atoms. In princi- 
ple, each Mn impurity atom contributes one hole, but 
heavy compensation intrinsically present in the system 
leads to a "dilute" carrier system with rii 3> rth, where 
Ui and rih are the active local moment density and car- 
rier density respectively. Throughout this paper, we re- 
fer to the carriers as "holes" and the magnetic local 
moments as "impurities" or "Mn atoms" without any 
loss of generality. The origin of this heavy compen- 
sation is essentially unknown although As antisite and 
Mn interstitial defects, both acting as electron double- 
donors, are generally thought to be responsible. There 
is wide consensus that the ferromagnetic ordering of the 
impurity moments at low enough temperatures is caused 
by the local exchange interaction between the impurity 
atoms and the holes. This exchange coupling produces 
a long-range effective ferromagnetic interaction between 
the impurity local moments, leading to the ferromagnetic 
state of the system at temperatures T < Tc, where Tc 
is the ferromagnetic transition temperature. The pre- 
cise mechanism for the carrier-mediated magnetic cou- 
pling between the impurity atoms is still being discussed 
and debated in the literature, particularly for the metal- 
lic DMS systems where both free-carrier weak-coupling 
RKKY-Zener interactiorki*^ and strong-coupling Zener- 
double-exchange interaction^ have been invoked depend- 
ing on whether the holes reside in the semiconductor 
valence band or in the Mn-induced impurity band re- 
spectively. The issue of whether the carriers reside in 
the semiconductor valence band or in the impurity band 
in the semiconductor band gap has been controversial.-^ 
The hole binding energy (~ 150 meV above the valence 
band edge) contributed by Mn in GaMnAs is fairly deep 
in the band-gap, and therefore impurity band physics is 
certainly operational in the strongly localized insulating 
DMS regime where the physical situation must be of car- 



riers strongly bound to randomly localized magnetic im- 
purities in the system. In our earlier publications^&iSiS 
we have developed a polaron percolation theory for DMS 
impurity ferromagnetism in the strongly localized insu- 
lating regime by showing that the localized holes produce 
magnetic clusters of bound magnetic polarons with the 
polaron size increasing with decreasing temperature un- 
til the magnetic clusters overlap through the whole sam- 
ple leading to the ferromagnetic temperature-induced po- 
laron percolation transition (more details on this scenario 
are provided in Sec. IIII|I . Such a magnetic polaron tran- 
sition scenario has later been explicitly verified in direct 
numerical simulations of the DMS Hamiltonian.^" In the 
current work we generalize our polaron percolation pic- 
ture by explicitly including in the theory the direct (i.e. 
not carrier-mediated) antiferromagnetic exchange inter- 
action between the impurity atoms which play an impor- 
tant role in II- VI materials, where the direct antiferro- 
magnetic exchange interaction is strong and at larger im- 
purity concentrations in III-V materials, when the mag- 
netic impurities are more likely to be at short distances 
from each other. 

The DMS materials are interesting systems from the 
fundamental perspective of frustration physics because 
there are independent competing ferromagnetic and anti- 
ferromagnetic contributions to the magnetic Hamiltonian 
arising from totally distinct physical processes: the ferro- 
magnetic contribution is carrier-mediated and effectively 
long-range whereas the antiferromagnetic contribution is 
extremely short-ranged direct Mn-Mn exchange. The dis- 
order inherent in random locations of the impurity atoms 
(i.e. Mn in Gai-a^Mn^jAs) coupled with the competing 
long-ranged (carrier-mediated) ferromagnetic and short- 
ranged direct antiferromagnetic interaction among the 
impurity atoms leads to interesting frustration physics 
in the system. This frustration will eventually lead to 
a non-ferromagnetic ground state (possibly, a spin glass 
state) in the system at high enough impurity concentra- 
tion and we study this DMS frustration physics within 
our polaron percolation model. We believe that this ran- 
dom competition between ferromagnetic and antiferro- 
magnetic coupling leads to the disappearance of DMS 
ferromagnetism at high impurity concentration and is re- 
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sponsible for the generic absence of ferromagnetism (with 
a few exceptions with very low values) in II- VI DMS 
materials. As described in the rest of this paper, we find 
that Tc shows a maximum as a function of the magnetic 
impurity density rii due to this random competition, and 
for large enough impurity density the system has a non- 
ferromagnetic ground state. We would like to mention 
in this context Ref. ITll in which a spin glass state was 
observed in an insualting (Ga, Mn) N alloy (which is a 
III-V dilute magnetic semiconductor). It is possible that 
the spin glass freezing is due to the random antiferro- 
magnetic coupling between magnetic impurities, which 
becomes dominant at high impurity concentration. 

The outline of the paper is as follows: The model is 
introduced in Sec. ^ In Sec. IIIII we describe the per- 
colation theory approach we use to deal with the prob- 
lem at hand. In Sec. IIVI we investigate suppression of 
the ferromagnetic transition temperature by direct an- 
tiferromagnetic interaction between the Mn impurities. 
The suppression is due to the fact that the coupling of 
spins to charge carriers and the strength of the inter- 
action between spins strongly depend on the concentra- 
tion of the magnetic atoms. We show that the approach 
used in Sec. IIVI fails at very large concentrations of Mn 
atoms. This case of large Mn concentrations is consid- 
ered in Sec. where a different technique is used. The 
results of the paper are summarized in Sec. IVII 



II. THE MODEL 



The magnetic Hamiltonian of the system has the form 



the case for the exchange coupling between a hole and a 
Mn moment in GaMnAs. 

The direct antiferromagnetic interaction between mag- 
netic impurities also decays exponentially with the dis- 
tance, 



J^^(r)=J„^Fexp(-9 



(3) 



and the corresponding decay length a is much smaller 
than that of the impurity-hole interaction, a <^ L, since 
the former is determined by the (small) decay length 
of the electron wave function of an inner shell of the 
magnetic impurity, while the latter is determined by the 
(larger) decay length of the wave function of an impurity- 
band carrier localized at some defect in the crystal struc- 
ture. 

Parameters of Hamiltonian (^), Jq, Jq^ , a, and L, to- 
gether with impurity and hole concentrations ni and rih 
make the full set of parameters determining the proper- 
ties of the system. In a realistic diluted magnetic semi- 
conductor system, these parameters obey the following 
relations: 

> — TTT > a , L > a . (4a) 



1/3 



1/3 



The fact that the charge carriers in the system under 
consideration are localized means that 
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The relations between the pairs of parameters 
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where indices j and k run over magnetic impurities and 
localized holes respectively, Sj and denote spins of 
magnetic impurities and localized holes, Jirjk) is the 
constant of impurity-hole exchange interaction, with its 
dependence on the distance rjk between the hole local- 
ization center and the magnetic impurities given by 



J(r) = Joexp ( -- 



(2) 



with L being the hole localization radius. The precise 
origin for the hole localization is irrelevant for our the- 
ory - it could be Coulombic binding to the impurity ac- 
ceptor or disorder-induced Anderson localization or any 
other relevant localization mechanism. The sign of the 
impurity-hole exchange interaction does not matter for 
the physical properties of the impurity subsystem, since 
the effective exchange interaction between impurities in- 
duced by the first term of Hamiltonian is always fer- 
romagnetic (see, for example, Ref. 5). We choose Jg > 
(that is antiferromagnetic) in our consideration without 
any loss of generality simply because that happens to be 



are not predetermined in our consideration. 

In this work, we will treat both impurity and hole spin 
dynamics classically. For impurity spins, it is justified by 
their relatively large value 5 = 5/2 in real experimental 
systems, such as GaMnAs. The hole spins, even though 
not so large (s = 3/2 in GaMnAs), are strongly coupled 
to the magnetic impurities closest to the hole localization 
centers, which makes them "massive" (that is effectively 
classical) as far as their interaction with more distant 
impurities is concernedi^ In any case, the quantum prop- 
erties are unlikely to be of qualitative importance in de- 
termining Tc. 

For analytic description of the system described by 
Eqs. CJ-Q), it turns out to be important whether re- 
lation 
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(6) 



is obeyed or not. The former case is more interesting 
physically, so it will be considered in this paper. The 
latter case requires a completely different formalism to 
treat but yields qualitatively similar results. In what 
follows we confine ourself to the former case, assuming 
that the inequality © holds. 
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III. PERCOLATION THEORY APPROACH 
A. Basic principles 

To understand the physics of ferromagnetic transition 
in a diluted magnetic semiconductor system, it is easier 
first to consider the system at temperatures well above 
the ferromagnetic transition temperature, and then fol- 
low the evolution of the system as the temperature is 
lowered. So at high enough temperatures (T ^ Jq) spins 
of holes and magnetic impurities are not correlated. As 
the temperature goes down, magnetic impurities close 
enough to hole localization centers tend to align their 
spins in the direction opposite to that of the hole they 
are close to. The resulting complex of a hole and mag- 
netic impurities polarized by it is called a "bound mag- 
netic polaron." The radius of a bound magnetic polaron 
grows as T — > 0. When two or more bound magnetic po- 
larons overlap, their spins must have the same direction. 
These overlapping polarons make polaron clusters; the 
ferromagnetic transition occurs when a macroscopic "in- 
finite cluster" appears. In our earlier works^ we con- 
sidered this ferromagnetic transition in the absence of 
antiferromagnetic interaction between magnetic impuri- 
ties (J(^^ = 0). While this antiferromagnetic interaction 
certainly suppresses ferromagnetic transition in the sys- 
tem (to what extent, we will see later), the basic me- 
chanics of the transition remains the same: bound mag- 
netic polarons grow, merge, make polaron clusters, and 
finally formation of an infinite cluster signals ferromag- 
netic transition. We note that the DMS magnetic polaron 
percolation transition is not the usual concentration- 
driven T = percolation transition, but a more subtle 
temperature-driven transition involving the coalescence 
of magnetic polaron clusters. 

Thus a bound magnetic polaron is the basic unit of the 
mechanism of ferromagnetic transition in diluted mag- 
netic semiconductors with localized charge carriers. We 
demonstrated in our earlier paper® that in the absence 
of antiferromagnetic interaction between the impurities 
they can be integrated out provided Eq. ^ holds, and 
the effective energy of the interaction between two po- 
larons with distance r between their centers can be ob- 
tained: 



27r 

fcff (?-,cos6',r)|jAF^o = ^ Y 



L^TOi In ( — 
Jo 



x:|exp(-^)cos0 (7) 

where 9 is the angle between the spins of two polarons. 
One can see that the interaction between the polarons de- 
cays exponentially with distance between them. This re- 
duces the problem of ferromagnetic transition in a system 
of bound magnetic polarons to the analogous problem of 
randomly placed spins with ferromagnetic exchange in- 
teraction between them decaying exponentially. 

The latter problem was considered by Korenblit et ali^ 
in the context of ferromagnetic phase transitions in di- 



luted magnetic alloys. They proposed mapping of the 
physical problem at hand to the mathematical problem 
of percolation transition of same-radius randomly-placed 
spheres. Such a mapping is possible because of two fac- 
tors: the exponentially fast decay of interaction between 
the spins and the characteristic distance between spins 
being much larger than the decay length of the interac- 
tion [see Eq. Ij4b|l ] . Under these conditions, the distri- 
bution of the nearest-neighbor couplings is wide. This, 
in turn, means, that an any given temperature T only 
a small fraction of nearest-neighbor-coupling strengths 
is of the order of T, and all the others are either much 
stronger or much weaker. The basic idea of the approach 
of Korenblit et alm^ is that we can postulate all spins with 
coupling strengths between them stronger that T to be 
locked in the same direction, and completely neglect all 
the other couplings, which, in general, are exponentially 
weaker. The justification for this approach is that there 
are very few couplings of the order of temperature, and 
their influence on the general picture is negligible; the 
validity of this justification will be discussed in the next 
subsection. 

The formal mapping to the problem of percolating 
spheres is done in the following way:^'^ all spins are re- 
placed with spheres of the same radius R, which is half 
the distance between two spins whose coupling exactly 
equals temperature. When two spheres (i.e spins) are 
closer than 2R, they overlap (and the coupling between 
the corresponding spins is larger than T). The spheres 
that overlap make clusters (which are magnetic clusters 
discussed above), and the ferromagnetic transition tem- 
perature is reached when the sphere radius (which grows 
as the temperature goes down) becomes large enough for 
the infinite cluster to appear. This problem of overlap- 
ping spheres has only one parameter, nR^, where n is the 
sphere concentration, and can be easily solved numeri- 
cally. Mapping the results back to the original physical 
problem of interacting spins, we get expressions for the 
ferromagnetic transition temperature Tc and other phys- 
ical characteristics such as magnetization, susceptibility 
etcj^iM In the case of the system described by Hamil- 
tonian ^ with J^^ — the transition temperature is 
given bji^ 



sS\Jo\ [Lnfl^^^Jm/nh 



exp 



0.86 



Ln 



1/3 



(8) 



We note that the solution of the problem involves two 
distinct mappings: first, the mapping of the DMS bound 
polaron Hamiltonian of holes and impurities to the prob- 
lem of random spins with ferromagnetic exchange, and 
then the random spin problem to the mathematical per- 
colation problem of random spheres. We also note that 
there is no characteristic temperature other than T^. in 
this scenario. 
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B. On the applicability of the percolation picture 

In this section we discuss the hmits of appUcabiUty of 
the percolation theory to the ferromagnetic transition in 
strongly disordered systems. Despite being an extremely 
useful tool in dealing with strongly disordered systems 
(sometimes, the only tool applicable), the percolation 
theory, as applied to temperature-driven phase transi- 
tions, has its drawbacks. Probably, the most important 
drawback, which is the center of this subsection's dis- 
cussion, is the inability of the percolation theory to ac- 
count for thermal fluctuations. We would like to em- 
phasize that this drawback is not specific to the perco- 
lation of bound magnetic polarons, which we consider 
in this and our earlier papers^*SiS In fact, any treatment 
of the temperature-driven ferromagnetic phase transition 
in a strongly disordered using the percolation theory, in- 
cluding the paper by Korenblit et al.^ in which this 
treatment was first proposed, inevitably neglects ther- 
mal fluctuations. In this section we argue that despite 
this problem, the ferromagnetic transition temperature, 
as predicted by the percolation theory, differs from the 
real transition temperature only by a numerical factor of 
the order of unity. 

We argue that the percolation picture of the 
temperature-driven ferromagnetic phase transition in a 
strongly disordered system of spins is qualitatively cor- 
rect. When the temperature goes down, ferromag- 
netic correlation is inevitably established first between 
the most-strongly-coupled spins, then spreading across 
weaker couplings. In such a system above the ferromag- 
netic transition temperature, one has ferromagnetically 
ordered regions (ferromagnetic clusters) without strong 
correlation between the magnetic moments of different 
clusters. As the temperature goes down, more couplings 
become saturated, so the new clusters appear and the ex- 
isting ones grow and merge. Clearly at some low enough 
temperature this growth and merging of clusters will re- 
sult in appearance of a correlated region which spans the 
whole sample (infinite cluster). This scenario is referred 
to, for example, in the context of the Griffiths phase 
physics, which was introduced before the first applica- 
tion of the percolation theory to the problem of ferro- 
magnetic phase transitions. It seems to be universally ac- 
cepted as far as the qualitative description of the physics 
is concerned. We have recently discussed the relevance 
of Griffiths phase to DMS systems in the context of the 
magnetic percolation transition in Ref. 0. 

However, it is hard to build a rigorous controllable 
approach based on the above scenario. A reasonable 
approximation^^ described in Sec. IIII Al is to assume 
that at any given temperature T two spins connected by 
some exchange coupling J are completely uncorrelated 
if J < T, and locked in the same direction is J > T. 
Clearly, thermal fiuctuations are completely left out of 
this approach. Despite this fact, such an approach should 
be adequate far from the transition point, when the clus- 
ters formed by the interacting spins have finite size. Any 



problems arising in estimating Tc can be taken care of 
by introducing numeric corrections of the order of unity 
to the physical quantities calculated. The infinite clus- 
ter is, however, a different matter. By definition, spins 
making the infinite cluster must all be correlated across 
arbitrarily large distances. Small but finite deviations 
from the strict alignment may accumulate over large dis- 
tances resulting in a complete loss of long-range coher- 
ence. Still one should agree that the exponent in Eq. (O 
is predicted by the percolation theory correctly. Indeed, 
there can be no ferromagnetic order at T 3> Tc with Tc 
given by Eq. (jS)) , because no connected network of cou- 
plings larger than such T would exist. The connected 
network of such couplings does appear at T ^ Tc. The 
transition temperature, however, is determined not only 
by the strengths of the couplings making the connected 
network, but also by the topology of the network itself. 
The latter, however, is characterized only by power (not 
exponential) dependence on the system parameters, as it 
follows from the scaling assumptions of the percolation 
theoryii^ Therefore the exponent in Eq. ||SJ) is the only 
exponent which may enter the expression for the true Tc- 
Now we are going to argue that not only the exponent, 
but also the parametric prefactor at the exponent for the 
true Tc is correctly predicted by the percolation theory. 
To begin with, let us note that the relation between the 

1 /3 

parameter p = of the percolation theory prob- 

lem and the temperature reads 

p = ny^ilog(^) , (9) 

with nl/'^L <C 1 [see Eq. (j4b|l ] . and the critical value of p 
being pc ~ 0.86. In the vicinity of the transition point, 
p Ki pc, we have 

b-Pcl ^ \T - Tc\ \T~Tc\ 

Pc ^ 0.86 Tc Tc ' ^ ' 

Therefore, at |T — Tc\/Tc ^ 1, that is away from the 
thermodynamic transition, we still have \p — Pc\/pc ^ 

1 /3 

L <C 1, so the system is still near the percolation 
threshold as far as the cluster topology and the scaling 
relations of the percolation theory are concerned. 

Now let us consider the system of spins at some tem- 
perature T = ATc, where A < 1 is a constant of the order 
of unity. Let us consider the spins making up the infinite 
cluster (in terms of the percolation theory) and argue 
that this infinite cluster must be in the ferromagnetic 
state at some value of A, which is not parametrically 
smaller than unity. As the basis for our arguments, we 
are going to use the "links, nodes and blobs" pictur o^^i^^ 
describing the topology of the infinite cluster in the vicin- 
ity of the percolation transition, see Fig.Q] According to 
this picture, the infinite cluster is a network of "links" 
connecting "nodes", with some occasional "blobs" em- 
bedded into the links. The link length is proportional to 
|p — Pc I to some power of the order of unity. One may ar- 
gue that the transition temperature of such a network is 
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do state that it is coherent with the established notions 
of the percolation theorjii^ while the contrary statement 
would conflict with these notions. The final answer may 
be given only by a numerical simulation, and we point 
out that the Monte-Carlo results-'^^ on the ferromagnetic 
transition temperature in diluted magnetic semiconduc- 
tors with localized charge carriers is in excellent agree- 
ment with our theory^, which is an indication that the 
reasoning of this section is correct. 




FIG. 1: Structure of the infinite cluster (after Refs. and 



of the order of Junk, where Junk is the effective exchange 
interaction between the nodes as it is transmitted by the 
link. The distribution of the couplings along the link 
is wide, the weakest coupling is exponentially small as 
compared to the average one, and the length of the link 
is large only as a power of the same parameter that en- 
ters the exponent. In such a system, it is the weakest 
coupling which gives the value of Junk, and therefore the 
critical temperature must be of the order of the weakest 
coupling in a link, which is of the order of as given by 
Eq. (jSJ. Therefore, the "real" critical temperature may 
not be parametrically smaller that the one predicted by 
the percolation theory. 

Now let us argue that the "real" critical temperature 
should not be parametrically larger than the percolation 
theory prediction. In principle, there is only one link 
of the strength of Junk between two nodes, so if it was 
the only link at all, the ferromagnetic temperature would 
have to be of the order of Junk- Even though other links 
connecting the nodes are weaker, they may significantly 
increase the critical temperature provided there is a suf- 
ficient number of them. We, however, argue that the 
number of weaker links is not large enough to compen- 
sate for their weakness. Indeed, the number of weaker 
links connecting two nodes may be as large as at most 
some power of the scaling parameter \p — pc\/pc, which 
at T ~ Tc is of the order of n^^i, see Eq. The link 

strength, on the other hand, is widely distributed due to 
the exponential decay of the interaction and the probabil- 
ity of parametrically many links having their strength of 
the order of Jimk is exponentially small, so the connec- 
tion between two nodes is dominated by the strongest 
link connecting them, whose strength is Junk: and the 
critical temperature may not be parametrically higher 
than Jiink- 

This reasoning leads us to the conclusion that while the 
ferromagnetic transition in the system under considera- 
tion is not, rigorously speaking, a percolation transition, 
the result for Tc given by the percolation theory may dif- 
fer from the real transition temperature only by a factor 
of the order of unity. We do not make any claims that 
our reasoning is a rigorous proof of our statement, but we 



IV. SUPPRESSION OF 

Because the decay length a of the inter-impurity an- 
tiferromagnetic interaction is much less than the char- 
acteristic distance nl^^ between them, see Eq. (|4a|) . the 
distribution of the coupling strength © between neigh- 
boring impurities is wide. Assuming the positions of the 
impurities to be random and uncorrelated, one can eas- 
ily calculate the probability P{J) that antiferromagnetic 
coupling J^^ between an impurity and its nearest neigh- 
bor is smaller than J: 



P{J) = exp 



3 



tAF 

In^^ 
J 



(11) 



It has been demonstrated'^'^'^^ that ferromagnetic in- 
teraction between two bound magnetic polarons (sepa- 
rated by the distance r) occurs mostly due to weakly- 
polarized lens-shaped region between the polaron cen- 
ters, with "thickness" L ln[(r/ Jq) exp(r/2L)] and "ra- 
dius" vTr, see Ref. In the expression ^ for the 
effective energy of the polaron interaction with J^^ — 0, 
the term in square brackets is essentially the number 
of magnetic impurities in this lens-shaped region, whose 
(weak) polarization by one polaron affects the other po- 
laron and vice versa. Non-zero antiferromagnetic inter- 
action suppresses this interaction in the following way: 
if the distance R between two impurities in the lens- 
shaped region is small enough, so that J^^{R) > T, 
the spins of these two impurities are locked in the op- 
posite directions, which means that they are not affected 
by the polarons and drop out of the group of impurities 
transferring interaction between the two polarons. The 
probability that an impurity is not excluded from this in- 
teraction transfer equals P{T) with function P given by 
Eq. (|ll|l . Therefore, the effective energy of interaction of 
two bound magnetic polarons in the presence of (direct) 
antiferromagnetic interaction between magnetic impuri- 
ties is given by Eq. Q with an extra factor of P{T) in 
the right-hand side, which accounts for the reduction of 
the number of impurities mediating interaction between 
the two polarons [factor in brackets in Eq. {Tj)]. 

Since the concentration of holes is relatively low, 
rihi^ ^ 1, the relative variations of £cS {r,cos6,T) for 
different pairs of neighboring polarons are very large. 
It enables us to use percolation theory to establish the 
temperature of ferromagnetic transition in the system. 
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as described in Sec. IIIII According to the percolation 
theory)^*i^ the ferromagnetic transition temperature 
is given by the condition 

feff (0.864, cos 0, Tc) =Te. (12) 

For not very large J^^\ Eq. ^ can be solved for Tc 
self-consistently by iterations, which yields: 

Tc ^ sS\Jq\ (Ln]l^\^/ni/nh exp ( - ^^/g ) 

/in ^ a^rij \ , , 

This expression for Tc becomes invalid when the antifer- 
romagnetic interaction between the impurities becomes 
so strong that most of them are locked that is when 
a^rii ^ Li^rih- The next section deals with this case. 



V. QUANTUM PHASE TRANSITION 

At some critical concentration of magnetic impurities, 
the Curie temperature is suppressed to zero by the di- 
rect antiferromagnetic interaction. This is a quantum 
phase transition point separating a ferromagnetic state 
and a disordered state in the system of magnetic impu- 
rities (with the density of magnetic impurities being the 
control parameter) . 

Let us briefly and qualitatively discuss the possible 
structure of the disordered state. It is reasonable to start 
with a model in which there are only randomly placed 
magnetic impurities with antiferromagnetic interaction 
© between them, and no holes. The classical ground 
state in the continuum model is not frustrated as the in- 
teraction between impurities is a well-defined (exponen- 
tial) function of the distance and the latter is a random 
variable. In the continuum model, the probability of find- 
ing, for example, a cluster of three impurities located in 
the corners of an equilateral triangle is exactly zero. One 
can say that the strong Possionian disorder in the contin- 
uum model lifts the degeneracy of the classical grounds 
state. 

At low temperatures, the quantum nature of the spins 
of magnetic impurities should reveal itself. To get an in- 
sight into the ground state of the system one can use the 
idea of the real-space renormalization group technique 
frequently used in the studies of low dimensional antifer- 
romagnetic systemSfi2: In this method, one systematically 
integrate out the strongest bonds gradually reducing the 
energy scale at each step of the decimation procedure. 
The strongest bonds are considered frozen into a singlet 
state. The locked impurity pairs are then removed from 
the system; however, the couplings between other spins 
get renormalized due to virtual triplet excitations across 
the strong bonds. The main assumption of the real-space 



renormalization group technique is that physical proper- 
ties are determined by an effective low-energy distribu- 
tion of couplings (which in some cases has a universal 
form essentially independent of the initial distribution). 
The averaged characteristics of the spin system (such as 
correlation functions) are determined by very rare events, 
corresponding to a fraction of spins locked into a singlet 
being separated by a very large distance. These rare 
bonds dominate long-range correlation functions. It is 
important to emphasize that typical properties of the 
system are drastically different from the averaged prop- 
erties and essentially determined by the majority of spins 
locked into a singlet with their nearest neighbours (one 
also should keep in mind that Griffiths effects are much 
less pronounced in higher dimensions). Therefore, the 
quantum ground state of our system can be represented 
as a set of randomly oriented spin pairs (valence bonds) 
locked into singlets. Such a state can be called a valence 
bond glass. Using this term, we can conclude that by 
increasing the density of magnetic impurities in a DMS 
system, one may induce a phase transition from a fer- 
romagnetic (low impurity concentration) into a valence 
bond glass state, in which the typical correlation function 
decays, at the lengths of the order of characteristic sep- 
aration rij between the impurities. At larger scales, 
the local impurity spin orientation can be considered as 
a random variable. 

The actual quantitative details for the quantum phase 
transition to the valence bond glass phase are beyond the 
scope of the current work, where we consider the (clas- 
sical) competition between random ferromagnetic and 
antiferromagnetic couplings induced by disorder at high 
magnetic impurity concentrations. It is clear, however, 
that the Tc-suppression given in Eq. (|13l) leads to a glassy 
magnetic state for large values of Ui, and quantum effects 
are obviously important in the elucidation of this glassy 
state. 



A. Bound magnetic polaron in a random medium 
of magnetic impurities 

Now let us consider one hole localized among magnetic 
impurities, still at zero temperature. Closer to the hole 
localization center, its interaction with the impurities is 
stronger than the antiferromagnetic interaction between 
the impurities, so all impurity spins in the vicinity of the 
hole localization center are aligned in the same direction. 
Away from the hole localization center, its interaction 
with the impurities becomes weaker, and at some dis- 
tance i?p the antiferromagnetic interaction between im- 
purities prevails over their interaction with the hole. To 
describe this behavior quantitatively, we introduce rif(r), 
which is the concentration of the impurities aligned by 
their interaction with the hole spin rather than by their 
antiferromagnetic interaction with neighboring impuri- 
ties. nf(r) equals the product of the impurity concentra- 
tion rij and the probability P[J(r)] that the coupling of 
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and impurity to its nearest neighbor is weaker than its 
couphng J(r) to the hole. Using Eqs. Q and l(TT)) . we 
arrive at 



rif (r) = 7ii cxp 



4 ^ / r\3 

-—nan 
3 



(14) 



The characteristic length i?p, at which 7if(r) decays is 
given by 



Rp — 



1/3 

n- a 



(15) 



d' 



and will be called the "zero-temperature radius of a 
bound magnetic polaron." 

The magnetic impurities surrounding a bound mag- 
netic polaron have some preferred spin directions deter- 
mined by their antiferromagnetic interaction with their 
neighbors. A bound magnetic polaron "feels" this as an 
effective magnetic field Brand, whose amplitude and di- 
rection are random and are determined by the (random) 
configuration of impurities around it. Now we will es- 
timate the characteristic magnitude of this field. There 
are two random contributions to this field. One, B|:^]j 
comes from the interaction of the hole with the mag- 
netic impurities which are not polarized by it. The other 

(2) 

contribution, 6];^^^^, is due to the interaction of the im- 
purities polarized by the hole with the unpolarized ones. 
The energy coming from interaction of an impurity spin 
polarized by the hole with other polarized impurities is 
also random, but it docs not depend on the hole's spin 
orientation, so it should not be included into the effective 
magnetic field which represents the action of the random 
impurity medium on the hole spin. 

Rigorous analytic evaluation of these two contributions 
is hardly possible. One could naively expect that they 
can be expressed in terms of integrals of n{{r), rii — nt{r), 
and J(r) over volume. However, while rif (r) is indeed the 
concentration of impurities whose spin direction is deter- 
mined by the hole spin, rij — nf(r) is not the concentra- 
tion of randomly oriented impurities that determine the 
magnitude and direction of Brand- The reason is that in 
addition to two above mentioned categories of impurities 
- the ones polarized by the hole and the ones belong- 
ing to the random medium which extends away from the 
holes - there is a third category. An example of impuri- 
ties of this third category is given by two impurities at 
the distance r^, < Rp from the hole localization center 
with the distance between them small enough for the an- 
tiferromagnetic interaction between them to be stronger 
than J(r*). One can easily see that these two impuri- 
ties do not belong to those making nf(r), but since this 
pair is able to rotate as a whole it will orient itself in the 
most energetically favorable position with respect to the 
spin of the hole, it will not contribute to Brand- We are 
not aware of any analytic way to separate these impuri- 
ties from the random continuum, and therefore we must 
limit ourselves to estimates of Brand instead of a rigorous 
evaluation. 



For this estimate, we simplify the picture to the fol- 
lowing: we assume that all the impurities whose distance 
to the hole localization center is smaller than Rp have 
their spins set by the hole spin, and that all the impuri- 
ties beyond Rp belong to the random continuum. In this 
"setup" the estimation for Brand is straightforward and 
yields: 



B. 



(1) 

rand 



IiLR'IJq exp 



Ru 



-°rand "i -"p-^O ^Xp 



1/3 



(16a) 
(16b) 



In these equations, the prefactors should probably be dis- 
regarded since they are likely to be artifacts of the sim- 
plification made above. The exponent Jq exp(— _Rp/L) = 
J(^^ exp(— 1/nV^a), however, which is notably the same 



for both B^^^^ and -Brand' highly likely correct, so it will 
be taken as our estimate for Brand. Since the method em- 
ployed above was rather crude, we may not be sure that 
a numerical factor of the order of unity [similar to 0.86 
in Eq. ((SJ] was not missed, so we write the equation for 
Brand the foUowing form: 



,(2) 



In 



B, 



and 



7AF 
•^0 



1 



1/3 

a 



(17) 



B. 



Interaction of two polarons in random 
continuum 



Now, having considered interaction of one bound mag- 
netic polaron with the random medium made by interact- 
ing magnetic impurities, we turn to interaction between 
two polarons surrounded by interacting magnetic impu- 
rities and separated by distance r from each other. In 
Ref. 5 and Sec. II VI of this paper we evaluated the inter- 
action strength between two bound magnetic polarons 
with high precision, that is up to a numerical prefactor. 
Such a precision would be excessive in the case considered 
in this section, since this interaction is to be compared 
with the interaction with the random medium surround- 
ing polarons, for which only estimate H17|l is available. 
Therefore we may limit ourselves to an estimate of the 
polaron interaction strength. The (ferromagnetic) inter- 
action comes from the interaction of each polaron's hole 
with the magnetic impurities polarized by the other, with 
the concentration of the latter given by Eq. H14I) . The 
strength of this interaction reads: 



£cs (r, cos 9, T) 



-Jo exp ( — - — I cos( 



(18) 



with all dimensionless prefactors omitted due to the rea- 
sons outlined above. Comparing Eq. H18|) with Eq. H17|) 
we immediately come to a (natural) conclusion that two 
polarons at zero temperature will have their spins aligned 
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FIG. 2: Temperature of ferromagnetic transition Tc as a 

1/3 

function of dimensionless impurity concentration nj a, with 
n^^L = const. The solid line is determined by Eq. 11311 . The 
zero-temperature transition point is found in Sec. [3 The 
exact behavior of the dotted curve are beyond the scope of 
this paper. 

in the same direction if the distance between them satis- 
fies the condition 

r < 2i?p , (19) 

where Rp is the zero-temperature polaron radius given 
by Eq. H15|l . and will have their spins pointing into two 
uncorrelated directions, which are determined by their 
interaction with the random medium of magnetic impu- 
rities if condition H19() does not hold. 

Employing notions of the percolation theory we used 
above, we can use parameter 



to characterize the ground state of the system. If 1 , 
the bound magnetic polarons form a ferromagnetically- 
ordered infinite cluster, and the ground state of the sys- 
tem is ferromagnetic. If 7^ <C 1, the interaction with 
the randomly-polarized medium of magnetic impurities 
breaks correlation between the cluster spins, and no long- 
range ordering of bound magnetic polarons is possible. 
The ground state of the system will be randomly ordered, 
with occasional isolated polarized regions around bound 
magnetic polarons and their finite-size clusters. Note 
that this statement is in agreement with result of 
Sec. lIVI which states that the ferromagnetic critical tem- 
perature becomes strongly suppressed (as compared to 



its value when J^^ = 0) exactly when V ~ 1. We delib- 
erately refrain from calling this random state a spin glass 
state since the words "spin glass" imply some very spe- 
cific properties, which this randomly ordered state may 
not possess. 



VI. CONCLUSION 

The results of this paper presented in Sec. IIVI and |Vl 
can be summarized in Fig. [3 At very low concentra- 
tion of magnetic impurities, the ferromagnetic transition 
temperature is low, because of the obvious reason that 
the impurity atoms mediate ferromagnetic interaction 
between bound magnetic polarons. Low concentration of 
impurities means weak interaction between polarons and, 
therefore, low transition temperature. As the concentra- 
tion of impurities goes up, the ferromagnetic interaction 
between bound magnetic polarons becomes stronger, and 
Tc increases. At lower concentration of magnetic impuri- 
ties, antiferromagnetic interaction between them can be 
neglected, due to its small decay length a [Eq. as 
compared to the large decay length L of the impurity- 
hole interaction. At a certain value of the concentra- 
tion of magnetic impurities, this antiferromagnetic inter- 
action becomes important and eventually the interaction 

1 /3 

between the polarons weakens, and Tc goes down as rij a 
grows, see Eq. ()13|l . We note that such a non-monotonic 
dependence of the ferromagnetic transition temperature 
was observed in experiment [see, e. g. Ref. At a cer- 
tain value of the concentration of magnetic impurities the 
Curie temperature is suppressed to zero and a quantum 
phase transition from a ferromagnetic state into a valence 
bond glass state is anticipated [see Sec. 0. The actual 
nature of this disordered valence bond phase remains an 
interesting topic for future theoretical studies. 

We have also argued in this paper that, although the 
DMS ferromagnetic transition is strictly speaking not 
a percolation transition, the polaron percolation theory 
should provide a reasonable estimate for the ferromag- 
netic transition temperature Tc, see Sec. IIIII 
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